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We analyse the out of equilibrium behavior of an Ising spin chain with an asymmetric kinetic 
constraint after a quench to a low temperature T. In the limit T — > 0, we provide an exact 
solution of the resulting coarsening process. The equilibration time exhibits a 'glassy' divergence 
: exp(const/T 2 ) (popular as an alternative to the Vogel-Fulcher law), while the average domain 



We show that the equilibration 



te 

length grows with a temperature dependent exponent, d ~ i 1 
time t eq also sets the timescale for the linear response of the system at low temperatures. 
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Even after decades of research, understanding the dy- 
namics of glasses remains a challenging problem (see 
e.g. One of the main features of glassy systems is 

that their relaxation time r increases quickly as the tem- 
perature T is lowered. A popular representation of this 
increase (for so-called 'fragile' glasses 0) is the Vogel- 
Fulcher (VF) law, r - exp[-const/(T - T )\. This pre- 
dicts that r diverges at temperature T , and the latter 
has therefore been associated with the temperature at 
which a true thermodynamic glass transition (achievable 
only in the limit of infinitely slow cooling) would take 
place. However, other functional forms for t(T) that 
have been proposed do not exhibit singularities at any fi- 
nite T, indicating the absence of a thermodynamic glass 
transition. Among these, the exponential inverse temper- 
ature squared (EITS) form r ~ exp(const/T 2 ) is popu- 
lar. Experimentally, it is difficult to distinguish between 
VF and EITS behavior due to obvious limitations on the 
longest accessible timescales; both can represent the ex- 
perimentally observed t(T) in many materials Q. Thus 
analytical results are desirable to shed light on this con- 
troversy. In this work we solve a simple dynamical model 
exhibiting glassy dynamics and find EITS behavior. 

To model relaxation in glassy systems theoretically, 
one can postulate some kind of quenched disorder, ei- 
ther in terms of some underlying microscopic Hamilto- 
nian (as is done in spin glasses) or more phenomeno- 
logically by making assumptions about the phase space 
of the system (e.g. in terms of hierarchical or ultramet- 
ric structures or energy barrier distributions 0,||). 
So far the main theoretical justification for either VF or 
EITS behavior comes from the latter approach; the EITS 
law, for example, is motivated by considering activated 
dynamics in a landscape of Gaussian distributed energy 
barriers 0. 

The alternative approach is to consider simple mod- 
els whose dynamics directly induce glassiness. Exam- 
ples include systems with kinetic constraints fi]] or en- 
tropic barriers ||, and driven diffusive models |10|. Such 
an approach is more obviously relevant to the dynamics 



of structural glasses (where quenched disorder is absent) 
since one does not need additional arguments that relate 
quenched and dynamically 'self-induced' disorder [ |TT| . 
The present work provides a first example where EITS 
behavior emerges directly from a microscopic model with- 
out imposed quenched disorder; instead energy barriers 
arise naturally from dynamical constraints. 

We consider a chain of spins in a uniform field, whose 
dynamics is nontrivial due to an asymmetric kinetic 
constraint. This model was introduced by Jackie and 
Eisinger [|l2| and has recently been rediscovered |13|. We 
study in particular the behavior after a quench to a low 
temperature T — * 0. We solve the resulting coarsening 
dynamics exactly in this limit and find two main results: 
Firstly, the equilibration time of the system diverges as 
t eq ~ exp(l/T 2 ln2) (EITS behaviour). Secondly, before 
equilibrium is reached, the average domain length grows 
as d ~ t T 2 , with an exponent that varies continuously 
with temperature. This novel anomalous coarsening is a 
consequence of the dynamical constraint, which produces 
scale-dependent energy barriers which grow as the loga- 
rithm of the domain size. Finally, we show that t eq is not 
just the timescale for equilibration after a quench, but in 
fact is also the timescale for relaxation of spin-spin cor- 
relations in equilibrium (at low T); this relaxation time 
therefore also has an EITS divergence at low T. 

The model comprises a chain of L spins Sj € {0, 1} 
where 1 < i < L; periodic boundary conditions imply 
that the left neighbor of si is s^. The dynamics for a 
given temperature T are defined as follows: At any time, 
only spins whose left neighbor is up (i.e., has the value 
1) can flip. For such 'mobile' spins, the rate for down- 
flips 1 — > is 1, while the rate for up-flips — » 1 is 
e = exp(— 1/T). Detailed balance is obeyed, and the sta- 
tionary distribution is the Boltzmann distribution for the 
trivial Hamiltonian H = 5^ i=1 Sj. For low temperatures 
the equilibrium concentration c = e/(l + e) of up-spins 
is small. Since these spins facilitate the dynamics, the 
system evolves slowly for small T. Moreover to eliminate 
an up-spin one first has to generate an adjacent up-spin. 
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Thus there are energy barriers in the system's evolution. 

We will be interested mainly in the behavior after a 
quench from equilibrium at some high initial temperature 
T\ ~ 1 to T <C 1. The basic objects that we use for 
the description of the system are domains. As shown 
by the vertical lines in . . . 1 1 000 1 1 1 1 1 1 01 1 001 1 1 1 1 1 01 1 . . . , 
a domain consists of an up-spin and all the down-spins 
that separate it from the nearest up-spin to the left. The 
length d of a domain also gives the distance between the 
up-spin at its right edge and the nearest up-spin to the 
left. Note that adjacent up-spins are counted as separate 
domains of length d = 1 . In equilibrium, the distribution 
of domain lengths and its average are 

P m (d) = e/(l + e) d , 4 q = l + l/e. (1) 

Now consider what happens after a deep quench to 
T<l,e< 1. The equilibrium concentration of up-spins 
at the final temperature T is c = 1 / d = e + C(e 2 ); hence 
the equilibrium probability of finding an up-spin within 
a chain segment of finite length d is 0(de) and tends to 
zero for e — » 0. In this limit (e — » at fixed d), the flip- 
ping down of up-spins therefore becomes irreversible to 
leading order. In terms of domains, this means that the 
coarsening dynamics of the system is one of coalescence 
of domains: an up-spin that flips down merges two neigh- 
bouring domains into one large domain. During such an 
irreversible coarsening process, no correlations between 
the lengths of neighboring domains can build up if there 
are none in the initial state Jl4| ]. For the present model 
the equilibrated initial state consists of domains inde- 
pendently distributed according to (0). Therefore a 'bag 
model' ]l4j or 'independent interval approximation' for 
the dynamics, which is defined by neglecting correlations 
between domains, becomes exact in the low-temperature 
limit (always taken at fixed d). 

We now estimate the typical rate T(d) at which do- 
mains of length d disappear by coalescing with their right 
neighbors. Because domain coalescence corresponds to 
the flipping down of up-spins, T(d) can also be defined 
as follows. Consider an open spin chain of length d, with 
a 'clamped' up-spin (sq = 1) added on the left. Start- 
ing from the state (so, si, . . . , sa) = 10 ... 01, T^ 1 (d) is 
the typical time needed to reach the empty state 10 ... 00 
where spin Sd has 'relaxed'. Any instance of this relax- 
ation process can be thought of as a path connecting the 
two states. Call the maximum number of 'excited' spins 
(up-spins except sq) encountered along a path its height 
h. One might think that the relaxation of spin s ( i needs 
to proceed via the state 11. . . 1, giving a path of height 
d. In fact, the minimal path height h(d) is much lower 
and given by 

h{d) =7i + l for 2"- 1 <d<2 n (2) 

where n = 0, 1, . . . This result is easily understood for 
d = 2™ jl5) . To relax the 2™-th spin S2» , one can first flip 



up s 2 n-i and use it as an 'anchor' for relaxing S2»- The 
corresponding path is (with s 2 n-i and S2" underlined) 
1...0...1 -> 1...1...1 -> 1...1...0 -> 1...0...0and 
reaches height h(2 n ) = h(2 n ~ 1 ) + 1; the +1 arises be- 
cause the anchor stays up while the spin 2™~ 1 to its 
right is relaxed. Continuing recursively, one arrives at 
h(2 n ) = h(l) + n; but h(l) = 1 because the only path 
for the relaxation of s± is 11 — > 10. To prove (||) more 
generally, define d(h) as the length of the largest single 
domain that can be relaxed by a path of height < h. 
Because of detailed balance, any relaxation path can be 
reversed, yielding a path of the same height from the 
empty state to the state 10 ... 01. In the same way, let 
us define 1(h) to be the maximal length of any spin con- 
figuration (ending in an up spin) that can be reached 
from the empty state by a path of height < h. One then 
has d(h+l) = 1(h) + 1 because to relax smh+i) one needs 
to flip up its left neighbor while exciting no more than h 
additional spins. A second relation is obtained from the 
relaxation of a configuration realizing the bound 1(h). 
Such a configuration contains h excited spins (due to its 
maximal length). To relax the first of these, no extra 
excitations are allowed (because of the ceiling h on path 
height); for the relaxation of the 2 nd , 3 rd . . . h-th spin, 
a maximum of 1, 2. . . h — 1 excitations are available. 
Summing the maximal length change at each step then 
gives 1(h) = ^2h~ = od(h' + 1). The above two recursions 
for 1(h) and d(h), combined with d(l) = 1(1) = 1, yield 
1(h) = 2 h - 1 and d(h) = 2' 1 " 1 , proving (§). 

At this stage we already see the key feature of the 
dynamics: the energy barrier for the relaxation of spin Sd 
is h(d) — l (the —1 comes from the one excited spin (sd) in 
the initial state). The rate for this relaxation is therefore 
T(d) = 0(exp[-(/i(d) - l)/T]) = 0(e h ^~ 1 ) g. Then 
eq. (||) tells us that the relaxation rate for domains of size 
d is T(d) ~ exp(— Ind/Tln2). Thus the energy barrier 
for the growth of domains increases logarithmically with 
domain size, giving a typical domain size growing as d ~ 
t Tln2 . Also, since d eq ~ exp(l/T) the equilibration time 
will grow according to an EITS law t cq ~ exp(l/T 2 In 2). 

From the scaling of T(d), the coarsening dynamics in 
the limit e — > naturally divides into stages distinguished 
by n — h(d) — 1 = 0, 1, . . . During stage n, the domains 
with lengths 2™ _1 < d < 2 n disappear; we call these the 
'active' domains. This process takes place on a timescale 
of 0(T^ 1 (d)) = 0(e~ n ); because the timescales for dif- 
ferent stages differ by factors of 1/e, we can treat them 
separately in the limit e — * 0. During stage n, the dis- 
tribution of inactive domains (d > 2 n ) changes only be- 
cause such domains can be created when smaller domains 
coalesce. Combining this with the (exact) bag model dis- 
cussed above, we have for d > 2" 

d T P(d,T)= P(d-d',T)[-d T P(d',T)} . (3) 

2"- 1 <rf'<2" 

The term in square brackets is the rate at which active 
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FIG. 1. Domain length distributions P n (d) at the end of 
stage n — 1 of the low T coarsening dynamics, for initial tem- 
perature Ti — oo. Open symbols and lines: Theoretical re- 
sults, calculated from (^), for n = (Oi initial condition), 
1 (□), 2 (O), 3 (A). Full symbols: Simulation results for a 
chain of length L = 2 15 and e = 10" 4 (n = 1, 2) and e = 10~ 3 
(n = 3). Inset: Scaled predictions 2 n_1 P n (d = 2 n_1 x) vs. a: 
for n = 1, . . . , 8. Bold line: Predicted scaling function (^). 

domains disappear; dl < 2 n because inactive domains do 
not disappear. We use the rescaled time r = te n ; during 
stage n of the dynamics and in the limit e— >0, it can take 
on any positive value r > 0. The initial condition for (j^) 
is the domain length distribution at the end of stage n—1 
of the dynamics, which we call P n (d) = P(d, r— >0). To 
calculate P n +i(d) = P(d, r— >oo), introduce the generat- 
ing function G(z, r) = Y^2 n - 1 <d T ) zd > an d its analog 
for the active domains, H(z,t) — 'Yl,2 n - 1 <d<2 n P{d 1 r)z d . 
From (§), one then finds 

d T [G{z, t) - H{z, r)] - -G(z, r) 9 T i? (z, r). 

This can be integrated to give [1— G(z, oo)]/[l— G(z, 0)] = 
exp[H(z,0) — H(z,oo)]. But at the end of stage n, all 
domains that were active during that stage have disap- 
peared, and so H(z, oo) = 0. Defining the initial con- 
dition for G as G n (z) = G(z,0) = T,2^<d P n( d ) z<1 
and similarly for the active generating function H n (z) = 
H(z,0), we then have finally 

G n+1 (z) - 1 = \G n {z) - 1] exp[H n (z)] . (4) 

This exact result relates the domain length distribu- 
tions P n (d) and P n +i(d) at the end of stages n — 1 and 
ii of the dynamics, as expressed through their gener- 
ating functions. Iterating it from a given initial dis- 
tribution Po(d) gives P n (d) for all n — 1,2,... We do 
this numerically by expressing (^) directly in terms of 
the probability distributions; the exponential is thus ex- 
panded into a series of convolutions of increasing order. 
Fig. [j] shows the results for the case where Po(d) is the 
equilibrium distribution (|l|) corresponding to an initial 
temperature of Ti = oo. Not unexpectedly, a scaling 



limit is approached for large n: The rescaled distribu- 
tions P n (x) = 2 n ~ 1 P n (d), where the scaled domain size 
is x = d/2 n ~ 1 , converge to a limiting distribution P(x) 
which is independent of the initial condition. Invariance 
under (^) gives an equation for the corresponding Laplace 
transforms g(s) and h(s) of P(x) 

g(2s) - 1 = [g(s) - 1] exp[h(s)} . 

We find a self-consistent solution 

00 I -nm-l r oo m , / m \ 

p(,) = E L 4r-/ n^ME— (5) 

m=l ' J1 r=l r \s=l / 

„ / si „ , s lrt(x — 1) 
= 0(a; - 1) - - Q(x - 2) — y - '- + ... 

X X 

where O(x) is the Heaviside step function. This series has 
singularities in the fc-th derivative at the integer values 
x = k + 1, k + 2, . . . The calculated P(x) agrees well with 
the previous results obtained by direct iteration of (^) 
(Fig. |l|). The average domain length in the scaling limit 
is given by d n = 2 n ~ 1 x; from the results for P{x) we find 
x = exp(7) = 1.78 . . ., where 7 is Euler's constant. 

In order to compare the results to simulations, con- 
sider starting from an equilibrated state at some initial 
temperature, say Tj = 00, quench the system to tempera- 
tures T <C 1 at time t = and observe its time evolution. 
If the results are plotted against the scaled time variable 
v = ln(<)/ ln(l/e) = Tint, then for T -> the n-th stage 
of the dynamics shrinks to the point v = n. In this limit 
we predict that, for n — 1 < v < n, the domain length 
distribution is P n (d) as defined by the recursion (Q). The 
average domain length d will follow a 'staircase' function, 
jumping at v — n from d n = ^2 d P n (d)d to d n+ \. In the 
large v scaling regime, this tells us that 2 u ~ 1 x < d < 2 v x 
(where x = 1.78 . . . from above), or \ < d/(xt Tln2 ) < 1 
when expressed in terms of ordinary time t. 

We can therefore say that the system coarsens with 
an exponent that depends on temperature and is given 
by Tin 2 to lowest order in T. By extrapolating this 
coarsening law to the equilibrium domain length d eq = 
exp(l/T) + O(l), we then also have that the dominant 
divergence of the equilibration time of the system for 
T -> is i oq = exp(l/T 2 ln2). 

In Fig. |^, we show the results of simulations for a range 
of values of e = exp(l/T). We used a waiting time Monte 
Carlo algorithm [p~7| combined with an efficient binary 
tree representation for the positions of the mobile spins. 
This let us access far larger systems (L = 2 15 ) and longer 
times (up to t = 10 10 ) than in previous simulations p3[ . 
The plateaus in d(v) that develop with decreasing e can 
clearly be seen, and their values are in good agreement 
with the predicted theoretical values. We also obtained 
the domain length distributions on the plateaus, by tak- 
ing data at the minima of (d/ dv)d{v) w.r.t. v. These are 
shown in Fig. ^ for the cases n — 1, 2, 3 and are again in 
good agreement with our theory. 
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FIG. 2. Evolution of average domain length d after quench 
from Ti = oo to T at t = 0, plotted on a log scale vs. v = T In t. 
Simulation results for four values of e = exp(l/T) are shown, 
obtained from a single run for a spin chain of length L — 2 15 . 
Bold line: Theoretical prediction for T — » 0. Inset: Theory 
for larger v and v — > oo asymptotes. 



Our result for the equilibration time t, 



cq — 

exp(l/T 2 In 2) is based on the extrapolation of the finitc- 
d coarsening behavior d ~ t T 2 into the equilibrium re- 
gion d = 0(l/e), where it is no longer strictly valid. We 
now show, however, that the same timescale is obtained 
from the initial decay of the spin-spin correlation function 
at equilibrium at low temperature T. It turns out that 
due to the asymmetric constraint the correlation function 
is site diagonal, ((sj(0) — c)(sj(t) — c)) = Sijc[R(t) — c] 
E3,E8[. Here R(t) is the probability that an up-spin at 
t = is also up at a later time t. With increasing t, 
it decays from R(0) = 1 to the equilibrium concentra- 
tion of up-spins, c = e/(l + e). To find the initial de- 
cay of R(t), consider again timescales t — <D(e~ v ) for 
finite v and e — ► 0. For v = n + 0, all domains of 
length d < 2" will have disappeared. Therefore only 
up-spins that bounded longer domains at t — will have 
an 0(1) probability of still being up. From the equi- 
librium distribution (|I|), one sees that they constitute a 
fraction (1 + e)~ 2 of the up-spins at t = 0, and hence 
R(y = n + 0) ~ 1 - 2 n e + 0(e 2 ) @. Neglecting cor- 
rections of 0(e 2 ), the quantity — lnR(v) thus lies be- 
tween 2"~ 1 e and 2 u e (for v > 0). Reverting to ordinary 
time, we have 1/2 < -[lnR(t)]/(t/t cq ) Tln2 < 1 for short 
times (t/t eq ) Tln2 <C 1. The relevant timescale that enters 
here is exactly the equilibration time t cq = exp(l/T 2 In 2) 
found above. We can thus identify the equilibration time 
for coarsening after a quench, with the equilibrium relax- 
ation time; both have an EITS-divergence at low T. 

Finally, we discuss briefly the spin-spin autocorrela- 
tion function for longer times (t/t eq ) Tln2 = 0(1), where 
the analysis becomes more involved |l2|,|l5|,|l8| . We have 
tackled this problem by extending the concept of domains 
to that of 'superdomains' which are bounded by up-spins 
that remain up on a given timescale. Combining this with 
a plausible hypothesis for the behavior of the relaxation 



timescales T~ 1 (d) for d = 0(l/e), the following scenario 
seems likely In the limit T — > 0, R{t) first decays 
linearly with the rescaled time variable 5 = (i/i eq ) Tln2 . 
This is compatible to lowest order with a stretched ex- 
ponential relaxation. But then the decay becomes much 
faster, and R actually decays to zero at a finite value of 
8. (For nonzero T, there is a crossover into a slower de- 
cay, presumably exponential in t, at late times.) It would 
also be of interest to study the relaxation times of similar 
models in dimension D > 1 |2lJ] . 
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